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The first passage time problem for Brownian motions flitting a 
barrier tias been extensiveiy studied in the iiterature. In particular, 
many incarnations of integral equations which fink the density of 
the hitting time to the equation for the barrier itself have appeared. 
Most interestingly, Peskir (2002b) demonstrates that a master inte- 
gral equation can be used to generate a countable number of new 
equations via differentiation or integration by parts. In this article, 
we generalize Peskir's results and provide a more powerful unifying 
framework for generating integral equations through a new class of 
martingales. We obtain a continuum of Volterra type integral equa- 
tions of the first kind and prove uniqueness for a subclass. Further- 
more, through the integral equations, we demonstrate how certain 
functional transforms of the boundary affect the density function. Fi- 
nally, we demonstrate a fundamental connection between the Volterra 
integral equations and a class of Fredholm integral equations. 



1. Introduction. Let {Wt)t>o be a standard Brownian motion started at zero and 
h : (0, oo) — > M be a continuous function satisfying 6(0) < 0. Define the first passage 
time (from above) of Wt to tlie curved boundary 6(t) to be: 

(1) T = M{t>0;Wt<b{t)} , 

with distribution function F{t) = P(t < t). The first passage time (FPT) problem seeks 
to determine F when b is given, while in the inverse problem we look for b given F. We 
will assume that b{t) is a regular boundary in the sense that P(r = 0) = 0. Sufficient 



conditions for regularity are given by Kolmogorov's test (see e.g. Ito and McKean (1965) 
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pp. 33-35). Furthermore, we allow 6(0) = — oo but we assume that whenever this is the 
case then there exists e > such that b is monotone increasing on (0, e]. 

The FPT problem for Brownian motion has a long history and available closed form re- 
sults appear to be sparse and fragmented. The few special cases include the linear bound- 
ary, quadratic boundary (see Salminen (1988) and Groeneboom (1989)| ), and square-root 
boundary (see Breiman (1967) Ricciardi, Sacerdote, and Sato (1984)| and Novikov, Fr- 
ishling, and Kordzakhia (1999)| among others). The celebrated method of images allows 
one to, at least theoretically, solve the problem for a class of boundaries, ba{t), which 
are solutions, for each fixed t, of implicit equations of the type 



-uba(t)~u'^t/2 



Q{du) 



where a > and Q is a positive cr-finite measure (see Lerche (1986) |. In this case the 
density function fait) is given by 



^(&a(t)M)/o°°^0(^)QW 
2t3/2 f^ ^i^'-^)Qi^de) 



fait) 



One way to tackle the FPT problem is to derive equations linking b and F. This is 
one of the primary motivations for studying integral equations in the context of the 



FPT. Peskir (2002b) presents a unifying approach to the integral equations of Volterra 
type arising from the FPT. Furthermore, the author generalizes the Volterra equations 
of the first kind. These equations are difficult to solve analytically but they are useful in 
a number of areas including the study of the small time behavior of F ( Peskir (2002a) [ ), 



numerical procedures yielding approximate evaluations of F (Park and Schuurmann 



(1976)1 [Smith (19727| [Durbin (197T) among others) or closed and semi-closed form 



approximations (e.g. Ferebee (1982) and Park and Paranjape (1974) |. Integral equations 
of Fredholm type are also useful in deriving known unique integral transforms of F (e.g. 
Shepp (1967)1 [Novikov (1981)[ ) and expansions of the FPT density (e.g. [Daniels (2000)| ). 

The first kind Volterra or Fredholm equations mentioned above can be viewed as being 
of the form EigiWr,T)) = g(0, 0) and as such are a direct result of the optional sampling 
theorem applied to an apropriate martingale giWs,s). Using this simple martingale 
result, our main aim is to present a unifying approach to the integral equations arising 
from the FPT and generalize the known class of integral equations. In Section 2 we 
examine such classes of martingales and provide a class of integral equations which 
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generalize all previously known Volterra integral equations of the first kind. Furthermore 
we examine necessary and sufficient conditions for the existence of a unique solution to 
a subclass of these equations. In addition, we outline a method, based on the method 
of images, for deriving new integral equations of Volterra type. In Section 3 we apply 
a similar martingale approach to derive Fredholm type equations. These equations are 
then used to provide concise alternative derivations of known closed form results for the 
linear, quadratic and square- root boundaries. Finally, we show the equivalence between 
the Fredholm and Volterra equations of the first kind for a particular class of boundaries. 

2. Volterra Integral Equations. The motivation behind connecting the martin- 
gale theory and the construction of integral equations for Brownian motion is perhaps 



best illustrated by the following well known Volterra equation (Peskir (2002b) I : 
(2) f\fy-b{s)\ F{ds) _ 1 y 







Vt — S J \/t — S \/t \y/t/ 
where (j) is the standard normal density function. The equality holds for all y < b{t) for 
continuous regular boundaries b. This equation can be written as 

(3) E{Xrl{T < t)) = Xo 

where the process X^ is defined as X^ = (j) (^^=f ) foi' fixed t > 0. Replacing (/> by the 
standard normal cdf, produces another well known equation which holds for all y < h{t) 
when h is continuous. Noting that Xg is a real- valued martingale for s < t and that 
Xtl{T > t)) = a.s., equation Q can be viewed as a product of the optional sampling 
theorem applied to the process XsAt and the stopping time r. 

Thus, the first step is to look for a class of martingales of the form Xs — m{Ws, s), s < 
t, satisfying E{\Xr\l{T < t)) = Jq \m{b{s), s)\F{ds) < oo and such that lims|( E{Xs1{t > 
s)) = 0. Suppose such a martingale exists and take a localizing sequence of stopping 
times s At, s < t. Then, applying the optional sampling theorem to X and s At and 
passing to the limit s | t, we obtain 

Xo = limEiXsAr) = limE{Xrl{T < s)) = f m{b{u),u)F{du) 

s^t sit Jq 

by the use of the almost sure identity Wr = b{T). The class of functions m for which the 
process Xg satisfies the above properties is a rather large class. A subclass of positive 



functions m can be constructed using the following classical result due to Widder (1944) 
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Theorem 1 ( Widder (1944) ) Let u be a continuous, non-negative function on I = 
(0,(5) X M, < 5 < oo. The following statements are equivalent: 

1) u satisfies the diffusion equation = | Uxx on I and lim(g .j:)^(o,e) u{s,x) = for all 
e < 

2) There exists a positive a-finite measure Q on [0, oo) such that u can be represented 
as 

(4) u{s,x) = ^0 Q{dB) . 

Given this result, define m(x,s) = u{t — s,y — x) for any t > and y < b{t). Tlien 
m satisfies tlie diffusion equation = —\mxx using the first part of the Theorem [l| 
Furthermore, the process Xs — m{Ws,s), s < t is a martingale (we can check directly, 
by computing the double integral, that S(|Xs|) = Xq < oo for all s < t). Checking the 
first condition, E{\Xt-\1{t < t)) < oo, we have: 

C\m{bis),s)\Fids)= f ^=^(^M_t^) F{ds)=u{z,t) 

Jo JQ \/t - S \ y/t - S J 

using equation ([2]). Furthermore, note that on the set {r > s} we have Ws > b{s). Take 
So close enough to t and such that for all so < s < t we have 6(s) > y. Such so exists 
since b is continuous and b{t) > y. Then 



E{Xs1{t > s)) = QidO) E [l(r > s)ct> ^ ) 



=u{t - s,y - 6(s)) . 

Taking the limit s ^ t and using the limitimg behavior of the function u as given in 
Theorem [T] above we see that limg^t E{Xs1{t > s)) = 0. This gives us the Volterra 
equation of the first kind: 

(5) u{t, y) = r u{t -s,y- b{s))F{ds) 

Jo 

for any y < b{t), t > 0. 

The integral representation of the function u (equation Q) is computable for several 
specific "degenerate" cases, such as when Q{d6) is a sum of Dirac measures or a uniform 
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measure over a compact domain. However, we have found one other general class of 
measures which lead to tractable forms for u itself, specifically when Q{d9) = 6~P~^d9 
for p < 0. In this case by direct calculation (see (|43|)) we have 

u{t,x;p) = e"''/'^'^'^Dp{-x/Vi)V~f^~\{p)/V2^ 



where Dp is the parabolic cylinder function (see Section B.l ). Note that even for p > 0, 
this particular u{t,x;p) still satisfies the diffusion equation ut = \uxx] furthermore, 
?x(t, x; 0) ~ (\){x I \fi) I \/t and u{t,x; —1) ~ ^{x/\/t) which are the kernels of the two well 
known Voltera equations. These observations motivate us to examine the function 

... ! ^ e-^-^Dp{{x-y)/Vt^s) 

(6) m(s,x;p) = {^ _ g^^(p+i)/2 ' P^V^x 

(for a fixed t > 0) more closeljQ We now proceed to show that m(s, Ws), s < t, is an 
honest martingale and derive a Volterra equation with kernel ni{s,b{s);p) by applying 
the optional sampling theorem. 

Define the process Xg = m{s,Ws;p), the stopping time Tt = t A t, and fix y G 
(— co,6(t)). We will use the optional sampling theorem on Xg, s < t, and Tf - it is 
important to point out that here time flows with s, while t represents a fixed time point. 

First we show that {Xs}s<:t is a martingale. Using the second order differential equa- 
tion (42), to which Dp is a solution, it is straightforward to show that 



(7) nis = -^rrixx ■ 

To check the integrability conditions, consider E [|m(s, VF5;p)|l(|Ws| > a)], a ^ y, s < t. 



Using the asymptotic behavior of the parabolic cylinder function (see ( 49 ) and ( 48 ) ) we 
obtain: 

E[\m{s,Ws;p)\l{\Ws\ > a)] 

\m{s,x;p)\ — dx + / \m{s,x;p)\ — dx 



-oo v27rs J a v2vrs 

^ ^dx+ -—^—^^dx<oo 



2tTS Ja {t - s)P+l/2^/2^ 



^We first came across this function through an alternative route prior to realizing the connection to 
the Widder's (1944) result. In fact, it is not apparent how Widder's result applies when p > 0. However, 
Dp {p > 0) can be written as a linear combination of parabolic cylinder functions with p < 0, the 
coefficients of which are space and time dependent. 
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Furthermore, m{s,x;p) is a continuous function in x on [—a, a]. Thus, 

E{\X,\) =E\m{s,Ws;p)\ 

=E[\m{s,Ws;p)\li\Ws\ > a)] + E[\m{s,Ws;p)\l{\Ws\ < a)] < oo. 

Therefore, Xg is a martingale for all p,y £ M. For s < t the process is a real valued 
martingale while for s > t it is a complex valued martingale. 

As before, on the set {r^ > s}, we have Wg > b{s) which implies (Wg — y)/\/t — s > 
{b{s) — y)/\/t — s — > oo as s t i because of the continuity of and the condition 
y < b{t). Thus, choosing sq close enough to t and such that for all sq < s < t we have 
b{s) > y and using the asymptotic behavior of the parabolic cylinder function (49), we 
obtain 

\Xs\l{Tt>s) = m{s,Ws;p)l{Tt> s) 

e \Dp{{Ws-y)/Vt^)\ 



= l(n>s) 
~ l(rt > s) 



(i_ s)(p+i)/2 

e ^(*--) {\Ws-y\/Vt^)P 
{t - s)(P+i)/2 



In particular, m{t, b{t);p) = lim^jt rn^s, b{t);p) = 0, for all p since b{t) > y. Furthermore, 
E{l{Tt> s)\Ws-y\n < / xP ^ dx 



b{s)~y y2'KS 



'l-KS Jb(s)-y 

Therefore, 

g 2(t-s) 

\\m.E\Xs\l{Tt > s)< Clim- = 

since — y) — > (K^) ~ y) > 0- Furthermore, whenever t > t, then IX^-J = \Xt\ = 
limg-ff = since t > t implies {Wg — y) > (6(s) — y) for all s < t and b{t) — y > 0. 
Therefore, 

E{\Xr,\)=E{\Xr\l{0<r<t)) 

= E{\m{T,b{Ty,p)\l{T<t))= I \m{s,b{sYp)\F{ds) . 

Jo 
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We already have that m(t, b{t);p) = for all p. Since m is a continuous function (since 
b is continuous) it follows that the last integral above is finite provided that 

\m{s,b{s);p)\F{ds) < oo 





for some small positive e. This is the case when 6(0) > — oo since m{0,b{0);p) < oo, 
so let us assume that 6(0) = — oo. Choosing e small enough and using the asymptotic 
behavior of the parabolic cylinder function we have 

(8) r\m{s,b{s);p)\F{ds) ~ P Ms) - y)-P-'\F{ds) 

Jo Jo 

and for p > —1 the last integral is finite since (6(s) — y)^^^^ ^ as s | and equals 



F{e) for p = —1. The case p < —1 follows from Lemma A.l and (41) in Appendix A. 
Therefore, for all p G M and y < b{t), we have E\Xrt \ < oo and by the optional sampling 
theorem Xq = E^Xrt); we have then proved the following result. 

Theorem 2 Let (Wt)t>o be a standard Brownian motion with Wq = 0. Let b : (0, oo) i— > 
M be a continuous function satisfying 6(0) < 0. Let t be the first-passage time ofW to 6, 
and let F denote its distribution function. Then for all p £M and y < b{t) the following 
system of integral equations is satisfied: 

(Q) e-'^D,{-y/Vi) _ f\-m^ Dpms)-y)/Vt^s) 

tiP+^)/^ ~ Jo (t_s){p+l)/2 ^^^^^ 

where F is the distribution of t. 

The set of integral equations (([9])) reduce to a set of well known integral equations 
when p = n, a non-negative integer, in which case (|9| becomes 

e- i Hn{-y/Vii) _ f g„((6(8) - !,)/v/2(I^) 



Here, Hn are the Hermite polynomials of degree n (see (44)). These equations were 



derived in e.g. Peskir (2002b) among others. In the next section we examine the limit 
y I 6(i) which allows the density and boundary to be tightly bound via the integral 
equations without the appearance of the arbitrary parameter y. Afterwards, we provide 
a richer class of examples. 
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2.1. Passage to the limit. The next step is to investigate what conditions on the 
boundary b are necessary to allow the limit y | b{t) in ^ to be taken. This limit is not 
straightforward for all values of the parameter p. To see this let us compute the limit as 
y I b{t) in equation Q with p = n = 1 assuming b{t) is continuously differentiable on 
(0,00) and 6(0) < 0. 



(h(s)-fc(t))^ 

First, in this case, for to > there exists some e > such that e < e 2(t^s) 



for all 



to < s < t since lim^-ft 
F{ds 



b'{t).0 = 0. Then we have 



to 



< 



(b(s)-b(t))-' 
2(t-s) 



to 



-F{ds) < 



2(t-s) 







(b('')-yr 



t g 2(t-s) 

liminf , F(ds) < liminf 

ym) \/t-s ym) Jo 



-F{ds) 

g 2(t-s) 

Vt-s 



-F{ds) 



< 00 



where the last equality follows from (|10[) with n = 0. Thus, when b{t) is differentiable 



t \b{t)-b(s)\ 

(i-s)3/2 



F{ds) < 00 since 



\b{t)-b{s)\ 
t-s 



is finite for all < s < 



/o ^ ^ ^^"^ therefore Jq 
t and in the neighborhood of the finiteness follows from Lemma A.l and (41). 

Second, for such boundaries, the corresponding density function of r is continuous i.e. 
F{ds) = f{s)ds where / is continuous on [0,oo) and /(O) = (see (Peskir 2002b) and 
dPeskir 2002"a| ). As a result. 



^3/2 



lim 

ytbit) Jo 



t (Hfh:y]l(b(s) - v) 



lim 

yW Jo 



{t - s)3/2 

ib(=)'yf (b(s) - b(t)) , , 

t-s 3/2 ym)Jo 



.(Hfkzy)! (b(t) - y) , , 



+ lim 2 

40 



l{u > z) exp ( -^u^ ( 1 + 



b{t - tz'^/u^) - b{t) 
z^ft 



f{t - tz^/u^)du 



where we have used the substitutions u 

U 17 1 ^ bU-tz^ /v?)-h(t) , 

above, l^or large u z, -^^ — ^ f 



b(t)-y 



and z 



^ in the third equality 



and thus there exists a positive constant 



o < 1 such that exp ^^^(l + 



< e ""^/2 for u ^ z. Therefore, since 
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/ is uniformly bounded, by the dominated convergence theorem we obtain 



Jo zio \ 2 y z\Jt 

rt (b(s)-b(t))'^ (h(s) - h(f)) I — 
since hm^jjo = 0. This last equality can be rewritten as 

(11) ^MMM = + r ^ (m-Ks)\ m -m 



The above equation was first derived by Ricciardi, Sacerdote, and Sato (1984) (see 
also Peskir (2002b) I . It demonstrates the complexity involved in exchanging the limit 
(as y I b{t)) and the integral in our new class of integral equations ^ - even for the 
"simple" case of p = n = 1. Nonetheless, we are able to compute this limiting case for 
a subclass of integral equations and the next result provides the required conditions on 
the boundary. 

Corrolary 1 Let (Wt)t>o be a standard Brownian motion with Wq = 0. Let b : (0, oo) 
M be a regular boundary and let r be the first-passage time of W below b, and let F denote 
its distribution function. Then, for all t > 0, the following system of integral equations 
is satisfied: 

c-TDp{-6(t)/i/i) _ f< DMH') - b(t))l./t^s) 



i) For all p < —1 when b is continuous on (0,oo) 

ii) For all —1 < p <0 when b is differentiable on (0, oo) 

iiijFor all < p < 1 when b is continuously differentiable on (0, oo) 

Proof. Note that Dp{x) > for all p <0. Define k{t) = hm^it ^''^/^^^ and g{s; t, y) = 

(M^)-y)^ 

e 4(t-s) Dp{{b{s) — y)/\/t — s). The function (7 is a continuous function in s on < s < t 

for all t > and y < b{t). Thus in order to apply the dominated convergence theorem we 

will show that g is dominated by an integrable function near s = and that g is finite at 
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s = t for all y < b{t). First note that when 6(0) is finite then \g{0; t, y) \ exists for all p and 
y < b{t) and when 6(0) = -oo then /q^^^ j^z^^TW^^id^) ~ /o^^^(Ks) " y)-''~^F{ds) < 
oo for some e(p) > and all p, y < b(t). The finiteness of the last integral follows 
from the fact that the integrand (6(s) — y)^^^^ is a monotone continuous function in y 
and thus for some y^ near b{t) it is dominated by (6(s) — y*)~^~^ which is integrable 
on (0,e(p)] by Lemma A.l . Thus we only need to show lims|t ^f}.lYp+i)/2 < cjo in or- 



der to apply the dominated convergence theorem since lims-ft ),'^lp+i)/2 — foi' V < K^)- 

i) Since lim,tt(t - s)~^p+^'^'^ = the case \k{t)\ < oo is straightforward. Suppose 

k{t) = oo. Then for s close to t, g{s;t,b{t)) ~ e~ '"'^t'-^)" ( ''^^)^*^ )p ^ using the 
asymptotic behavior of Dp{x) for large x (see ([49|)). Similarly, suppose /c(t) = — oo 
then the asymptotic behavior of c/(s; t, 6(t))(t - s)-(p+^)/2 ^(-^^^(^ _ g)-(p+i)/2 ^ 

(-^^^^^)-i-P(t - s)~(P+i)/2 ^ _ 5(s))-i-P I as s T i since - 1 > and 6 
is continuous. Therefore, taking the limit y | b{t) in ([9]), by the dominated convergence 
theorem the result follows. 



ii) We showed that when 6 is differentiable (and thus continuous) then //^ "7^^ 



< 



00, to > 0. Furthermore, differentiability implies k{t) = 0. Similarly as in part i) we 
see that Yun.s]tg{s]t,y)/{t — sy/'^ = for all y < 6(t) and and thus g{s;t,y)/{t — sy/'^ 
is bounded on [toi^]- By the dominated convergence theorem we can exhange the limit 



and the integral in 12 



iii) When 6 is continuously differentiable on (0, 00) then / is continuous on (0, t] for all 



t > (Peskir 2002b) and so Jq ^^^_fyj}+i)/2 d-s < 00 since < (p+ l)/2 < 1. Furthermore, 
\g{s; t, y) \ is bounded for < s < t for all y < b{t) since k{t) = 0. The result follows by 
the dominated convergence theorem. □ 

Note that the differentiability condition on the boundary in part ii) can be relaxed 
to \k(t)\ < 00 for all t > 0. In this case we still have //^ ^J^^ ^ ^0 > 0, using 
the same argument as before and the proof of part ii) is still valid. Also, it would be 



straightforward to extend the class of equations (|9j) and (12) to the class of equations 
with a complex valued parameter p. 
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2.2. Special Cases. For different values of p tlie parabolic cylinder function, Dp, can 
be written in terms of other special functions. The case when p is a negative integer 



covers the system of equations derived in (Peskir 2002b I as we will see in Case 3 below. 



Furthermore, equations (19]) and ( 12 1 can written in terms of the Whittaker function (see 



(45)) or confluent hypergeometric functions using the representation of the parabolic 



cylinder function for all values of p. When p is non-negative integer we already saw the 
connection with the Hermite polynomials which can be written in terms of the Laguerre 
polynomials. For p = —1/2 there is also a connection with the modified Bessel function 
of the third kind Kiy (see Case 4 below) . In this section we explore some of these special 
cases. 



Case 1: p = 



In this case ( 12 ) becomes 



g 2(t-s) 
^/t — S 



-F{ds) 



which can be written as 



(13) 



This equation was derived in (Durbin 1971 1 who uses a previous result by Fortet(1943). 



(Durbin 1971) uses the equation to obtain a numerical solution by approximating the 



boundary by straight line segments on subintervals (s, s + ds) and using available re- 



sults for crossing probabilities for linear boundaries. Subsequently (Smith 1972) recog- 



nizes ( |13| ) as a Generalized Abel equation and proposes Abel's linear transformation 
T : 5 — > g{t) I — tdt, to deal with the singularity of the kernel at (s = t). He then 
solves the equation numerically using standard quadrature methods. 



Case 2: p = -1 



In this case ( 12 ) becomes 
(14) 



Equation (14) was used in (Park and Schuurmann 1976) as a basis for numerical com- 



putation of the unknown density / using the idea of Volterra to discretize the equation 
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and solve the resulting system. This equation is especially attractive for numerical com- 
putations of / when b is given since the kernel K{t,s) = — b{s)) / \/t — s) is 
nonsingular in the sence that it is bounded for all < s < t. When h{t) = c < then 



(14) reads P{t < t) = 2<I>(c/\/t) which is the reflection principle for Brownian motion. 



Case 3: p = — n, n = 1, 2, 3... In this case \12\ becomes 

m - h{t)) 



-(b(t)-b(s))^ 
ft g 4(t-a) 

Vo'^ \~7W^s 



{t-s) — F{ds)= D^n 



We claim that (15) is equivalent to Peskir's system of equations (Peskir 2002b). Con- 



Gn{x). Using (51) we have 



sider the kernel of the integral equat ion {|15|) , wliicli is of the foriii — ^ ^^D—ji(^ — 



d_ 
dx 



Gn+l{x) = Gnix) 



and thus 



/X 
Gn{u)du + C 
-oo 



Taking x = and using (|52|) we see that C = 0. Therefore we can rewrite (|15|) as 
(16) 



Gn it-s)(--'y'F{ds) = Gnm/v~¥''-'^/' 

where n = 0,1,2,... and G„ satisfies the recursion formula Gn+i{x) = J^^Gn{u)du 
with Go{x) = <f>{x) since Do{—x) = e~^^/^. Therefore, the system of integral equations 



(15) is equivalent to the system of equations (16) which was derived in Peskir (2002b) 



This completes the proof of the above claim. 

The next two cases provide two new integral equations arising as specific cases of our 
general class: 



Case 4: p = —1/2. In this case, using (46), (12) becomes 



(17) 



1 



b(t) - b(s) _ (m-Hs))^ f (b(t) - 6 s M 
^ t-s 4t-s i ^ ^ 
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Case 5: A new class of equations can be derived from ( 12 ) using the recursive relation 
property (50) of the parabolic cylinder function. Using this relation and the class (12) 
for p < —1, we obtain: 



e 4t 



(b(s)-b{t))^ 

e 4(t-s) 

(t - s)P/2 

g 4(t-s) 
{t-s)- ^^A: 



V Vt-S / y/t-S \ \/t-S J) 



'(t_s)(p+2)/2 P- 



{b(s)~b(t))^ 
4(t-s) 



e 4t 



D,^,{-b{t)/Vt) + ^-^D,{-h{t)/Vt) 



{b{s)-b(t)y 



(t_s)(p+l)/2 \,/t^ 



-D 



P+i 



b{s) - b{t) 



+ 



Thus, from the last equality, we derive the class of equations for p < —1: 



(1 



{b{s)-b{t))^ 
4,(t-s) 



{t-s)(p+^y^ \ vi 



Ks) - bjt) 







For example, in the case p = —2 and using the results for p = and p = —\ together 

(b{t)-b{sy 



with (47), (18) becomes 
s 



(19) 







m - Ks) 



+ b{s)^ 



F{ds) = 



2.3. Uniqueness of a solution. Next we examine sufficient conditions for the bound- 



ary b such that the class of integral equations (12), which we denote as {i?p}p<i, has 
a unique continuous solution. We will first investigate the uniqueness of this system of 
equations for the case —1 < p < 1 and then generalize to the case p < —1. Suppose 
that b is continuously differentiable on (0,T] and assume lim^jo < oo for some 

e < 1/2. Note that lim^jo < cxd implies — oo < 6(0) < (since 6 is a regu- 

lar boundary and e < 1/2) and therefore the hitting density /(O) = (Peskir 2002a). 
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Therefore F{ds) = f{s)ds where / is continuous on [0, oo), . Denote (p + l)/2 = A so 
that < A < 1. Let 

52A+i(t) = e-'-^D2x+i{-b{t)/Vt)/t^ 

(b(s)-b(t)f 

K2x+i{t,s) = e D2x+ims)-Kt))/Vt^s) 
Then the class {-B2A+i}o<a<i of integral equations becomes: 

^* K2X+l{t, s) 

lo {t - sY 

Equations of this kind are also known as generalized Abel equations of the first kind. 
We know that the above equation has a continuous solution / ( Peskir (2002b) | ). Thus, 



(20) g.x^,it) = i '-^^^^f{s)ds 



to show uniqueness, it suffices to show that (20) is reducible to a Volterra equation of 
the second kind which has a unique solution. 



Using Lemma (A. 2) we have, |-Cp(^'^%=^^)| < Mp for some Mp > and all < s < 



t, p S R. Applying Abel's transform to equation (20) we obtain: 



(21) / = r r 



(U Jo J, (i.-()l-^(t-«)* 

where we have used Fubini's theorem (since |i^'2A+i| is bounded) to exchange the order 
of integration. Let gx{u) denote the left side of the above equation and K2x+i{u, s) = 
Is g)A dt = Jq '^^'^(Y-^)"-^^^^'^'' '^y then the last equation can be written as: 

(22) ~gx{u)= r K2X+i{u,s)f{s)ds 

Jo 



Next, we apply the standard technique of differentiation on u to reduce (22 ) to a Volterra 
equation of the second kind. First we show that 52A+i(0) = and that gx{u) has a 
continuous derivative for all ti > 0. For any Ai,A2 G M, since b{t)/\/t | — oo when t | 
for a regular boundary b, we have (using the asymptotic expansion of the parabolic 
cylinder function) 



/I X 2A1+A2 + I/2 



= lime- — H-b{t)r^-+'^-j =0 
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since oo > —b{0) > 0. In particular 52A+i(0) = 0. Also, using (51), 



(23) 



d92X+i{t)/dt = + Vtg2x+2{t) [^-^ - j 



^ _ Afi(2A+iW , 92X+2{t) , _ h{t)g2\+2{t) 
^ ' t ^ ' 2t 

Under our assumption on the boundary and since 6(0) > —oo each term in the last line 
goes to as t I and we obtain 

lim5'2A+i(*) = 0. 

Therefore, since h is continuously differentiable, it follows that 52A+i(^) ^^'^ 92\+i{'t) are 



continuous functions for all t > and since g2A+i(0) = 0, by Theorem 3, p. 5, Bocher 
(1909)[ gx{u) has a continuous derivative, for all u> 0, given by 



Next we compute the derivative, w.r.t. u, of the righthand side of (22 ). Since |I?2A+i((6(s) 
b{t))/\/t — s)\ < Cx for some Cx > and all < s < t it follows that i^2A+i(?/(^ ~ s) + 
s, s) < Ca for all < s < n and < y < 1, while Cxj^z^r^dy = CxB{l - A, A) ~ 
here B{-, •) represents the Beta function. Thus, by the dominated convergence theorem,: 

K2X+i{u,u) = lim^2A+i(^i,s) 



/•I 1 

/ hm K2X+1 {y{u - s) + s, s) J- rj- 

Jo 4m (1 — y) 

^2A+i(0)5(l-A,A) 



for all > 0. 
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Furthermore, using Lemma (A.2|, 
y^^'^-'\dK2x+iiyiu -s) + s, s)/du\ 



(b(s)~h(y{u-s) + s)) 



4y{u — a) 



p. ( b{s)-b{y{u~s)+s) ^ 

^''^'^ Vviu-s) \ y{u-s) 
{u — s)'^+^/2 y{u — s) + s 



X {b'{y{u — s) + s){y{u — s) + sY 



{b{y{u -s) + s)- h{s)){y{u - s) + sf 
2y{u - s) 



\Hl{y,u,s)\ 



< 



M 



for some constant M > and for all < s < u and < y < 1. Also 

M 1 



(M-s)^+l/2(l-y)l-AyA-(l/2-e) 

= MB(l-A,A-(l/2-e)) T 

Jo 



dy ) f{s)ds 

m 



-ds < oo 



{U - s)l/2+€ 

since / is continuous on [0,T]. Thus the derivative dK2x+i{y{u — s) + s, s)/du is domi- 
nated by an integrable function. Denote 

Kix+iiy^'^^s) = K2\+i{y{u - s) + s + 5y,s) - K2\+i{y{u - s) + s,s) 

for (5 > and note that /CgA+iC^' ^' '^)/'^ ~^ dK2\+i{y{u — s) + s,s)/du uniformly on 
(s, y) E [0, u] X [0, 1] as I 0. Therefore, if /i denotes the measure on [0, u] with Radon- 
Nykodim derivative / and denotes the measure on [0, 1] with Radon-Nykodim deriva- 
tive 1 ^(1 — y)^^-'*y'*'^(i/2-'=) ^ by Fubini's (applied twice) and the dominated convergence 
theorems, we have 



ru rl 

lim / / 

HO Jo Jo 



S (1 _ yy-XyX-{l/2-e) 



dyds 



-5^0 J[o,„]x[0,l] ^ 
[0,u]x[0,l] HO 6 

H\iy,u,s) , 

, .+1/2 <^ ^ ^) 

[0,«]x[0,l] [u - s)^+V^ 

Hl{y.u,s) 



-d{fj, X u) 



-d{n X u) 



u / ri 



(1 -y)l-AyA-(l/2-.) ^u-sY+y 



lis) 



zds 
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where the quantity 

is bounded by MB{1 — A, A — (1/2 — e)) and continuous for all < s < u with a possible 
discontinuity at u = s. Therefore the derivative of (22) w.r.t. u is given by 

ru+S 



~i f ^ y I So K2x+i{u + 5,s)f{s)ds - K2x+i{u,s)f{s)ds 
g\{u) =lim — 

]^ ru+S _ 

= lim- / K2\+i{u + 5,s)f{s)ds 

+ hm r K2X+,{u + 5,s) -K2X+i{n,s) ^^^^^^ 



<5iO Jo 

Thus we obtained the Volterra equation of the second kind: 

(25) -Ji^^ ^ /(..) + /" ^ai^iL^/(.).„ 

K2X+i[u,u) Jo K2x+i{u,u)[u- sy+^/^ 

Since R\ is finite on < s < u with a possible discontinuity only along the curve s = u 
and since / is continuous on [0, oo), by Theorem 3, p. 19, Bocher (1909) (25) has a 
unique continuous solution. Thus, the continuous solution to (20) is unique. Therefore, 
for each —1 < p < 1 equation Bp has a unique continuous solution. Furthermore, suppose 
p < —1 and g : [0,t] M is any continuous solution to Bp. Then, using the integral 
representation of Dp (and the fact that Dp{x) > 0, \/p < 0) and Fubini's theorem, we 
can write equation Bp as 



e V ^ 



Jo 

2 

OO pt { b{s)-h{t) ' 



\/t — s Jo \Ji 







nf^ f b{s) — h{t) . u \ /o ( \ POO — P — 1 

yt — s JO vt 



t _( Hs)-b{t) ^^^\ ,2 -{u/Vi~b(t)f/{2t) 

e V ^/ ' ds = ^ 

^/t- S xjt 

where the last equality holds for all n > and follows from the uniqueness of the Mellin 
transform. Taking the limit u | in the last equality it follows that any continuous 
solution to {i?p}p<_i is also a solution to Bq which has a unique continuous solution. 
Thus we proved the following result: 
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Theorem 3 For each T > let b(t) be a regular boundary, continuously dijferentiable 
on (0,T], and satisfy \b'{t)\ = 0{t~^) for some < e < 1/2. Then r, the first-passage 
time of the standard Brownian motion Wt to b{t), has a continuous density function, f, 
given as the unique continuous solution of the class of integral equations 

e-'-^D,{-b{t)/Vt) _ D,ms) - b{t))/Vt^s) 

tiP+^y^ ~ Jo {t-s)(P+^)/^ ' 

where p < 1. 

2.4. Functional Transforms. Next we consider some functional transforms of the 
boundary and the corresponding density functions. The new density functions can be 
easily expressed in terms of the original boundary and its density function using equation 
Bq and Theorem |3j Suppose b satisfies the hypotheses of Theorem 3 with corresponding 
density function / and introduce the functional transforms: 

(26) {Ti.b){t) = b{t) + at, a G M 

(27) {Tlb){t)=b{^t)/^, 7>0 

(28) {T,^.bm = (1 + [3t)b (y^) ' /3 > 

Note that we can set /3 < with t < —1/(3 in the last transform. Moreover, (Ti.6), (T2.6) 
and {T^.b) all satisfy the hypotheses of Theoremjsj Denote with /i, /2 and /s, respectively 
the corresponding density functions of the first-passage times of Wt to these boundaries. 
Using equation Bq we can easily find the relations between / and fi, i = 1,2,3. 
For /i we have: 



-{b{t)+aty/{2t) ft g 2(t37) 



\/t Jo \/t- s 



yft Jo Vt 



fl{s)ds 



s 



Therefore, due to uniqueness of solutions (Theorem |3]), we must have 

(29) fi{t) = f(t)e~^m~^'t/2 ^ ^-a{T-'\m) f^t) 

This result can alternatively be obtained by a simple measure change argument. 
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For /2 we obtain: 

-Hit) /(^■yt) ft p 2-t(t-s) 

f2{s)ds 



A/t Jo \/t - s 



U Jo \/U — X 

and, applying Theorem |3] 

(30) f,(t) = ^f{jt)=j'/\Tlf){t) 

This result can alternatively be derived through the time change t — > at. 
For /3 we obtain: 



Vt 



U Jo \Ju — X 



— g 2(l-;9ti) / g 2(u-i:) g l-/3u ^2 1{\-l}x)l (ix 

"v/m 7o ^Ju — X 

-b{uY/{2u) ^ 



u Jo V' 



'U — X 



0b'' (x) / T \ 



where we have made the substitutions u = j^^, x = jq^g^- Therefore, by Theorem jsj 
(31) /3(t) = / (^) exp (-/3(1 + /3t)62(^)/2) (1 + /?t)"3/2 



This result was obtained by Alili and Patie (2005) for more general boundaries, using 
probabilistic arguments. 

It is instructive and pleasing that the integral equations lead to a unifying derivation 
of all of these transformation results. Combining the three transforms into the single 
transform 

(T.6)(t) 4 ^-^^b [ytJ^) + « e /3 > 0, 7 > 
we have the following result for the corresponding density function fx- 
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Theorem 4 For each T > let b(t) be a regular boundary, continuously dijferentiable 
on (0,T], and satisfy \b'{t)\ = 0{t^'^) for some < e < 1/2. Let f be the continuous 
density function of the first-passage time of the standard Brownian motion Wt to h{t). 
Let b{t) = {T.b){t). Then the first-passage time to the boundary b{t) has a continuous 
density f given by: 

(32) f{t) = ^-^^TO*^ ^-{l+P^t)b{^){f3b{^)/2+a/^)-aH/2 

(1 + /37t)3/2 

Proof. Since b{t) is continuously differentiable then so is b{t) and thus the first-passage 
time of Wt to b(t) has a continuous density function. Application of the transforms 



, T2 , Tf succesively to b{t), transforms f{t) to f{t) and (32) follows from Theorem 

El □ 

3. Fredholm Equations. Similarly to Section 2, in this section we examine the 
well known martingale g-^'^s-^^Vs -^vhich gives rise to a Fredholm integral equation of 
the first kind. This equation is used to obtain alternative derivation of known closed 
form results for the linear, quadratic and square-root boundaries. Furthermore, as we 
will see, this equation is simply the Laplace transform of the integral equations ([9| 
for a particular class of boundaries. We assume that b is continuous on [0,oo). Let 
Ta = inf {t > 0; Wt < b(t) -\- at} with cumulative distribution function and define 
the set 

■^b{t) — {a G M; ba{t) = h{t) + at>c for some c < and all t > 0} 

Under the measure P* given by P*[A) = fj^Z{uj)dP{uj) where Z = \-aWt ^ 

distribution F. Then the equality Ep*{lr^<t) = Ep{lr^<tZ) becomes 

F{t) = f' Ep{e~^+"^'\Ta = s)Fa{ds) = rSp(e"^--e-'^^+°(^'-^=Va = s)Fa{ds) 
Jo JO 

*e"(^W+-)e-"^+^F„(ds) = t e^'^^^+'^Faids) 



Jo 

where we have used the almost sure equality Wt^ = b{Ta) + ar^- Since the above is true 
for all t > we have 

(33) Fa{dt) = F((ii)e-"^W-'^ 
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Under the assumption ba{t) > c , we know that < a.s. for all a £ A, where 
t'^ = mi{t > 0]Wt < c}. Since t'^ is almost surely finite then so is r^, which implies 



Fq(oo) = 1 and thus, for a £ Ah(^t)i using (33), we obtain the Fredholm integral equation 
of the first kind 



(34) 

with kernel K{a, s) 
and a G 



F{ds) = 1 



-Qb{s)- 



. Note that equation (33) holds for any boundary h 
while equation (34) holds for any a E -4.^(4). The latter equation can also be 



derived using the martingale property of the Geometric Brownian motion together with 



the optional sampling theorem and has been found as early as Shepp (1967) 



Next we extend equation (34) for complex values of a. Consider the processes 



Xt 
Yt 



-xWt~ 



* cos{yWt + xyt) 

* sin{yWt + xyt) 



for X, y G R. Both processes are martingales for all real x, y and thus the process Zt = 
Xt — iYt = g-'^^t-—t jg ^ complex valued martingale where a = x + iy. Define the class 
of continuous functions h. 

B = {h{t)] h{t) + ut > c, some c < 0, Vu < 0, for large t} 

Note that b(t) S B implies b{t) is uniformly bounded below and thus the corresponding 
first passage time is almost surely finite. 

Theorem 5 IfbGB and is continuous on [0, oo), then for all complex a with \ arg(a)| < 
7r/2, the equality 



(35) 



holds. 



'F{ds) 



1 



Proof. First notice that for b £ B equation ( 34 ) holds for all real a since b(t) + at > 
b(t) — \a\t > c for t large enough. We first look at the quantity £'(e^"^/^), r > 0. For any 
such r, since b G B, there exists an < N{r) < oo and a c G M such that for t > N{r) 
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we have ^/rh{t) > rt + = \fr{c + \/rt). Then 

£;(ei^) = ^(ei^l(r < N{r))) + E{ei''l{T > N{r))) 

< ^(ei^l(r < N{r))) + S(e-"^+v^^(^)~5^1(T > N{r))) 

poo 
/•oo 

Jo 

Next we apply the optional sampling theorem by showing £^(1^,-!) < oo and lim^^oo E{Zt'ir>t) 
0. For rc > and using the finitenes of -E(e'"^/^), r > 0, we have 



where c' is the uniform lower bound of b. Similarly, for x > 0, we can find an M{x, y) such 

2 2 

that for t > M(x, y) we have the inequality xb{t) — ^ ~^ t > 0. Then, for t > M(x, y) we 

2_ 2 

have \Xt\l{T > t) < e-^*W+^^*l(r > t) < l(r > t) and thus Wmt^^ E{\Xt\lr>t) < 
lim^^oo > t) = since r is almost surely finite. For x = 

E {\Xt\l{T >t))<E (ey''/H{T >t)) <E (^ey'/^n{T > t)) = j"^ ey''/^F{ds) < oo 

and thus liuit^^ E {\Xt\l{T > t)) < limt^^ Jl^ ey^'/'^F{ds) = 0. Thus, for all x > 0, by 
the optional sampling theorem, Xt and r satisfy E{Xt-) = Xq = 1. The same arguments 
apphed to the process Yt yield E{Yr) = Yo = 0. Thus E{Zr) = E{Xr) - iE{Yr) = 1 and 
the proof is completed. □ 



The above result gives an extension of the Fredholm equation (34) for boundaries 
belonging to the class B. When — f < arg(a) < ^ it is sufficient that b{t) is uniformly 



bounded below for equation ( 35 ) to hold since y^ — x'^ < 



Fredholm equations of the first kind are notoriously difficult to solve (even in the case 
when there is a unique solution). The two general cases in which explicit results are 
available are equations with kernels of the form K{at) or K{a — t). In the first case we 
can obtain the Mellin transform of the solution and in the second the Laplace transform. 
Next we examine boundaries which give rise to such kernels. The following results are 



well known, however, here we demonstrate that they all follow from equations (34) and 



(35) and illustrate their importance. 
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Example 1. 



b{t) = —a + bt, a > 0,b > 0. Thus A = {a > —b} and equation (34 1 becomes 

POO 2 

Jo 

If / is the Laplace transform of / then the above equation reads f{u) = e-b+\/iP+2u 
this is the Laplace transform of the well known Bachelier-Levy formula: 



^ -{a-btf/(2t) 



Example 2. 



b{t) = p^/i — q, q > 0, p / 0. Then A = {a > 0} and equation (34) becomes 

/■OO 2 

JO 

Multiplying both sides of the above equation by a^~^, x > and integrating a on [0, oo) 
we obtain 

roo /"OO 2 roc 

JO JQ Jo 

=> r f{t)dt r a^-'e-'^P^'-'^'da = ^ 
Jo Jo g^' 

"/(t)(2|)-tr(x)eP'/^Z)_.(p)dt= 

2 q^ 

t~-2fit)dt 



where D is the parabolic cylinder function. The last equality gives us the Mellin trans- 
form of / if we replace x with 2(1 — x), x < 1. Alternatively, by making the substitution 
i = e" in the last equation we obtain 



e-(f-i)«/(e«)j^ 



e 4 g ^ 



which gives us /(x), the Laplace transform of /(e"), after replacing x with 2x + 2 

-(2a;+2) 

(36) f{x) = --, X > -1 
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A similar approach was used in Shepp (1967) for the first-passage time to the double 



boundary ±a\/t + b. Novikov (1981) generalizes (36) for stable processes with a certain 
Laplace transform using a martingale approach. 
Example 3. 

b{t) = ^-q, p,q > 0. Take a such that 3f?(a) > 0. Denote a' = a{2pY/^ then 5R(a) > 
as well. Using a' equation ( 35 ) becomes 



'2 2 * 



f{s)ds 



-a q 



and after completing the cube under the integral, multiplying both sides of the equation 
by /3 > and integrating a over any contour C(a) with end points ooe~*3 and 
ooe*3 and | arg(a)| < 7r/3 (see Figure [I]), we obtain 

'^'K'b JC 



1 

27ri JC 

coo 



e-PP^e—f{t) 



2TTi Jc 



dadt 



c 



c 



The right hand side of the last equation is Ai{{2pY^'^[j3 — q)), where Ai is the Airy 



function (see (53)). Next we examine the contour integral on the left side. Define the 



contour C = C + ,r,^^i/3 and let zi , Z2 be points on C and their corresponding images on 
C be z[ and ^2 (see Figure Since the function under the contour integral is analytic, 
its integral over the simple closed contour ziz'iz'2Z2 is 0. Thus, sending zi to ooe's and 
Z2 to ooe 3 we obtain 



„/3(2p)i/3„^_i„3^^ = Ai(/3(2p)V3) 



ic Jc 
since the contributions on the legs ziz[ and Z2Z2 diminish in the limit. Therefore the 
Laplace transform of e e f{t) is given by 



(37) 



-at 



e 6 



f{t)dt 



7v 21/3 N 



The last example is an alternative derivation of (37) which was first obtained by Salmi- 



nen (1988)1 using measure change, and independently by [Groeneboom (1989) | using a 
factorization of the density f{t), involving a Bessel bridge and a killed Brownian mo- 
tion. 
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3.1. Connections with Volterra Integral Equations. Finally, we discuss the connection 
between the Volterra integral equations of Section 2 (see ([9])) and the Fredhom integral 



equations studied in the previous section (see (34)) for a certain class of boundaries. 



Let b{t) be continuous and uniformly bounded below, i.e. there exists a constant c < 



such that for all t > 0, b{t) > c. Such boundaries satisfy (34 ) for all a > 0. Set y < c < 
and a = ^/2P, /? > 0. Then, multiplying both sides of (34) by ^2P+i/2/3Pe^\/^, p < 
we obtain the equation: 

(38) r e-^^0F2P+i/2/3Pe-^(^(^)-^)F((is) = ^2P+^/'^pPeyV^ 







For any —p,z,l3 > 0, we have, (see Gradshteyn and Ryzhik (2000) 7.728), 

roo g-/3a;^-{p+l) 
Jo 0F2P+1/2 



-e 4,x D2p+l{z/^/x)dx = P^e 



and therefore the integral equation (38) can be written as: 

(b(s)-j;)2 , , , 



-pt_ 



-lit 



(t-s) 



P+i 



dtF{ds) 



tp+1 



dt 



'p+i[ 



/t-s 



{t - s)P+i 



F{ds) 



dt 



-I3t 



(y) 



tP+i 



dt 



imsart-aap ver. 2008/08/29 file: FPT_paper.tex date: February 24, 2009 



26 



JAIMUNGAL, KREININ AND VALOV 



Substituting 2p + 1 with p, and invoking the uniqueness of Laplace transforms, allows 
us to identify the terms in square braces and results in the class of integral equations 



for boundaries b{t) > c, t >0. 

4. Conclusions. In the first part of this article, we developed a new class of Volterra 
integral equations for the distribution of the first passage time(FPT) of a standard 
Brownian motion to a regular boundary. This new class generalizes and unifies the class 
of all such previously known integral equations. Interestingly, this class arises through 
the optional stopping theorem applied to an interesting and new class of martingales 
generated by the parabolic cylinder functions. Through an Abel transformation, we were 
able to prove uniqueness of a solution to the integral equations. Based on uniqueness, we 
were then able to consolidate the derivation of the FPT distribution to a transformed 
boundary. 

In the second part of this article, we generalized a class of Fredholm integral equations 
to the complex domain. These equations were then shown to provide a unified approach 
for computing the FPT distribution for linear, square root and quadratic boundaries. 
We believe that the method can be more widely applied by searching for specific factor- 
izations of the kernel that produce known transforms such as Mellin, Laplace, Hilbert 
and so on. Finally we demonstrated that there is a fundamental connection between the 
Volterra and the Fredholm integral equations studied in this work. 

There are several directions remaining open for future research. 

• The first is clear but difficult: how can this larger (uncountably infinite) class of 
Volterra integral equations be used to extract the FPT distribution? One way is 
to explore the flexibility of the parabolic cylinder function and its connection to 
other special functions. Furthermore, the continuum of Volterra equations provides 
more flexibility for manipulation such as integration and differentiation w.r.t. the 
parameter p. 

• The search for new Volterra equations of the first kind is related to identifying 
analytical solution to the heat equation. The search for such solutions, which gen- 
erate kernel functions with known properties, is another topic for future research. 



(39) 




for alio < s < t, p < 1, y < c. This is precisely the class of equations Q for p < 1, y ^ c. 
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Any linear combination of solutions to the heat equation is also a solution and 
possibly, in the limit, one can obtain Volterra equations with more informative 
kernels. 

We saw that taking the limit y t h{t) in ([9|, with p = 1, produced the Volterra 
equation of the second kind (111. This motivates the investigation of this limit for 
the equations with p > 1. We suspect that in the computation of this limit for 
p > 1, we can obtain new Volterra equations of the second kind and such equations 
are known to exhibit unique solutions and are generally easier to deal with than 
the Volterra equations of the first kind. However, such equations would hold for a 
restricted class of boundary functions. 

The class of Volterra equations is also a useful tool for the inverse first passage 
time probem. Though, in this context, the equations are highly non-linear the 
generalization of this class provides flexibility for their manipulation which could 
extract new information. 



APPENDIX A: USEFUL BOUNDING LEMMAS 

In this appendix we provide two useful Lemmas which provide important bounds on 
the boundary and the density. 

Lemma A.l Suppose 6 : (0,T] — > M is an increasing continuous function on (0,e] for 
some < e < 1 with &(0) = — oo. Let h : —>■ M and h{x) = 0(e"^^) for large x > 
and some < a < 1/2. Define the first passage time r = {s > 0; Wg < b{s)}. Then 

(40) r \h{-b{s))\F{ds) < oo 

Jo 

where F is the distribution of t. 

Proof. Without loss of generality we can assume —h{t) ^ for t < e. Define the first- 
passage time rb(s) = {t > 0;Wt < b{s)} for fixed s < e. Since b{t) < b{s) for t < s < e 
then F{s) < Fr^^^^{s) = 2<^>{-b{s)/ ^/s) for all s < e. Let g{s) = h{-b{s)), s < e and 
fix si < e and 6 > he such that ka^sisi) = e"''^^*^^ — 6 > 0. Define Sn such that 
go62(s„) _ /j^^^(sj^) -)- fi§_ Since e"''^^*^ is monotone decreasing on (0,e) with e"''^^'^-' = oo 
then Sn i is a monotone decreasing sequence. Let 

oo 

95{s) = 5 ^ l(s < sn-i) + {ka,s{si) + 6)l{s < e) 

n=2 
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Then < \g{s)\ < Me"^'(*) < Mgs{s), s G (0, e], for some M > 0, and by the dommated 

convergence theorem and the definition of g, there exists an e > such that 

/■<= °° 
\gis)\F{ds) <M / gsF{ds) = V F(s„_i) + Me'^''^''^ F{e) 

^0 n=2 

<2Ms ^{b{sn)/v^) + c< ^^-^ J2 msn)/v^) + C 

1 V ^ 1 

n=l ^ 71=1 



TT ^ 

n=l 

<oo 



TT 



n=l 



where C = Me"'^^^'^'^^ F{e). The second inequality on the second line holds since ^(x) < 
4>{x), X > while the last inequality follows from a < 1/2. This completes the proof. □ 

In particular, for y G M, > 0, we have 



(41) 



{y - 6(s))^F(ds) < OO 



Lemma A. 2 Let b{t) he a continuously differentiahle function on (0,r] with -co < 
6(0) < satisfying lim^jo l^'(OI^'^ < c« for some < e < 1/2. Then, for all < s < t,: 



m - b{s) 



< C 



and 



m - b{s) 

t - s 



f <K 



for some positive constants C and K. 

Proof. Since h is continuous on [0, T] the results hold for all < s < t < T. Since b is 
differentiahle on (0,r] the results hold on the curve T > s = i > 0. We only need to 
check the case ,5 = ^ = 0. For s = and t | we have that 



lim 



b{t) - 6(0) 


= lim 


m - 6(0) 




m 


t 



e.l/2-e 



lim|6'(t)|ft^/2-. 

40 







and 



lim 

no 



m - b{o) 



t 



f = lim|6'(t)|t' < OO . 

no 



This completes the proof. □ 

APPENDIX B: SOME SPECIAL FUNCTIONS 

In this Appendix we collect some important and useful results for the parabolic cylin- 
der and Airy functions. For additional properties of the parabolic cylinder and Airy 
functions see [Erdelyi (1954)] and [CFadshteyn and Ryzhik (2000)[ 
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B.l. The Parabolic Cylinder Function. The parabolic cylinder functions are 
solutions to the differential equation 

(42) <g+J,+ i/2-f)„ = 0. 

The complete set of solutions are u{z) = Dp{z), Dp{—z), iz). 
They admit the following integral representation for p < 0: 



(43) Dp{z) = / e-^^-^ /^x-P-^dx . 

i [-p) Jo 

Furthermore, for special sets of parameters p, they reduce to other more well known 
special functions: 

z 



(44) Dr^iz) = 2-"/2e-^'/4i?„ 

(45) Dpiz) = 2'/^^P/'Wy,^p/,,_,/J'-]z-'/' 



z^ 



(46) D_,/^iz) = yW2^i 



2' 



1/4 



4 



(47) D_2{z) = e'"/'^ (e-''/'^ -V2^z^-z)) . 

Here, iJ„ is the Hermite polynomial of degree n, W is the Whittaker function, K is the 
modified Bessel function of the third kind. 

For large argument z, they admit the following asymptotic expansions: 



(49) Dp{z) ~ z^'e"^"/^, \z\ T oo, |arg(z)| < Sir/A . 



(48) Dp{z) ~ ^^e'^'Pz-P-^e''/^ \z\ T oo, 7r/4 < | arg(z)| < 57r/4 , 



They are closed under derivative and integral operations in the following sense: 

(50) ^^f(^) = ^zDp{z) - Dp+,{z) , 

(51) ^^e-^'I^Dp{z) = e-^'/^Dp+,{z) , 

(52) I e-^ ''D_p{z)dz = r(p/2 + l) = -^-(^+1)^ " 
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B.2. The Airy function. The Airy function (for complext argument x) has the 
integral representation: 

(53) Ai{x) = I e^'/'^-'-^'dt 

27n Jc 

where the integral is over a path C with end points ooe~'^*/'^ and ooe'^*/^ (see Figure 
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